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Abstract 

Markopoulou and Smolin have argued that the low energy limit of LQG may suffer from a conflict between locality, as defined by 
the connectivity of spin networks, and an averaged notion of locality that emerges at low energy from a superposition of spin network 
states. This raises the issue of how much non-locality, relative to the coarse grained metric, can be tolerated in the spin network graphs 
that contribute to the ground state. To address this question we have been studying statistical mechanical systems on lattices decorated 
randomly with non-local links. These turn out to be related to a class of recently studied systems called small world networks. We 
show, in the case of the 2D Ising model, that one major effect of non-local links is to raise the Curie temperature. We report also on 
measurements of the spin-spin correlation functions in this model and show, for the first time, the impact of not only the amount of 
non-local links but also of their configuration on correlation functions. 

PACS numbers: 04.60.Pp, 89.75.Hc, 05.50.+q, 64.60.Ht 


In loop quantum gravity (LQG), states are described 
by spin networks, from which non-local (NL) links can 
emerge due to the mismatch of micro-locality and macro¬ 
locality 01 H@j. In other words, two nodes which are lo¬ 
cal microscopically in the fundamental graph can turn out 
to be macroscopically far away from each other and con¬ 
nected by a NL link relative to the coarse grained metric. 
Ends of the NL links carrying gauge field representations 
may be interpreted as gauge particles, while the ends of 
those carrying no gauge fields can be some unknown neu¬ 
tral particles and hence serve as a type of source of dark 
matter J2lID]. It is found that the MOND force law may 
also be recovered on a regular Euclidean lattice randomly 
decorated with some NL links characterized by some distri¬ 
bution function P(r), which is the probability of having a 
NL link between two nodes separated by a distance r mea¬ 
sured by the regular lattice metricyj]. Since NL links corre¬ 
late the fields at randomly chosen far away points without 
violating the local conservation of averaged energy, they 
may explain the non-locality in quantum physics and more¬ 
over, the origin of quantum physics according to Nelson’s 
formalism of stochastic quantum theory fill |@]. 


Simulation of 2D Ising Spin System with Random NL links 

We adopt Metropolis Monte Carlo method and apply pe¬ 
riodic boundary condition. We mainly simulate a 20 x 20 
2D square ferromagnetic Ising model randomly decorated 
with NL links based on the uniform distribution probabil¬ 
ity. Note that rewiring is not allowed and any two spins 
cannot be directly connected by more than one link. We 
take the interaction coupling, J = — 1, and Boltzman con¬ 
stant, k a = 1, such that all the quantities presented are 
dimensionless. The spin-spin correlation function g (r) 
needs some clarification. We study the pair correlation of 
the system as a function of r, defined as the minimum num¬ 
ber of links connecting one spin to another. Because non¬ 
local links break the symmetry of the lattice, g (r) should 
be obtained, for each r, by doing both ensemble average 
and the average over the number of pairs of spins with the 
same distance r. Mathematically, this is described by 

g( r ) = X [( a * CT f)r - (to) <0j»r] /#pairs(*.j)r> 

where the () denotes ensemble average, the subscript r de¬ 
notes the distance between two spins in the pair, and # 
means ’’the number of’. 


To address these questions, we study the statistical be¬ 
havior of 2D Ising models, with randomly added NL links. 
Such a system is related to a class of recently studied 
systems—small world networks (SWN), which are very 
useful in describing many interesting social systems, such 
as railway network and disease spreading yUa]01 • How¬ 
ever, our approach and focus in this letter are different from 
those for the SWN in that 1) we add a given number of NL 
links randomly to the regular lattice with uniform distri¬ 
bution probability, and 2) we study and find the effects of 
not only the amount of NL links but also the configuration 
of NL links on the system’s behavior, especially on corre¬ 
lation functions, which are tightly related to the system’s 
dynamics. To do so, we use Monte Carlo simulation. 


Effects of the Amount of NL links 

Since we assume a uniform distribution probability for 
the NL links, the only control we have on NL links is the 
number of them. Therefore, it is natural to look at the ef¬ 
fects of NL links on the system’s behavior by varying their 
number. 

Fig. Uli compares five curves, which correspond to the 
specific heat of a 20 x 20 regular lattice with respectively 
0, 10, 20, 50 and 100 NL links. One can see the increase 
of the critical temperature with the number of NL links. 
This is so because, for a finite system, adding NL links ef¬ 
fectively increases the dimensionality of the lattice and of 
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FIG. 1: Gy-curves of respectively a regular 20 x 20 Ising spin 
system, the same system with 10, 20, 50, and 100 non-local links. 
The data of this plot is from Dr. Lixin Zhan with Wang-Landau 
method. 



FIG. 2: Plot of T c /N versus N, the system size, in the case of 
total connectedness. 


course the average valence of the nodes. This argument can 
also be verified by noticing the decrease of the peak height 
as the number of NL links increases. Due to the finite scal¬ 
ing effect, the peak rises with the system size, L, defined 
by L D = N, where N is the number of nodes of the sys¬ 
tem and D is the dimension. Now that D increases with 
the amount of NL links, but with N fixed, L must decrease 
and so does the peak height. The last effect shown in the 
plot is that the half width of the Cy -curve also increases 
with the amount of NL links. This is understood by not¬ 
ing that 1) the system with NL links has more states and 2) 
there are more non-frozen states in the region of T <T C as 
T c increases with the amount of NL links. Similar results 
are reported for SWNs on Ising model by Cai et. al. 01, by 
HerreroQ, and by Hastings[10 ] l . 

One may ask if there is an upper limit of the critical 
temperature when we keep adding NL links. Yes, there 


1 Hastings proposed for a different model, an approximated scaling relation 
between the increase of T c and the probability p of having NL links (long- 
range links in his work). His model can be comparable to ours when p —> 
0. We have not found a close match with his prediction, probably because 
we do not evolve sufficently large lattices to obtain such a scaling behavior. 



FIG. 3: Comparison between the correlation function of the reg¬ 
ular 20 x 20 Ising model and that of the same lattice with 20 
non-local links randomly distributed. 


is a maximum T c , which is the very critical temperature 
in the extreme case when the system is totally connected. 
In the stage of total connectedness, T c raises with the sys¬ 
tem size, which is now equal to N. This is illustrated in 
Fig. Q. where T c /N is asymptotic to unity in the limit 
of JV —* oo. One can locate T c = 364 for a totally con¬ 
nected system of 400 spins (having 79800 links), which 
corresponds to a 20 x 20 regular Ising system with 79000 
NL links. This also implies that T c , after a very short linear 
range, increases slower and slower to approach the maxi¬ 
mally allowed T c . 

It is more interesting to see how the system’s correla¬ 
tion function behaves in the presence of NL links. Fig. 0 
shows three pairs of correlation functions, each of which 
contains one for a regular lattice and one for the same lat¬ 
tice with 20 NL links. The upper pair compares the two 
at their critical temperatures, one at T c = 2.3 and one at 
T' c = 2.5 due to NL links. One can see that both g (r)’s be¬ 
have similarly and indicate correlation lengths of the order 
of L, although with different decay rates. By getting rid of 
the boundary effect, we have g (r) oc r~° ,2188 for the reg¬ 
ular lattice, and g (r) oc r °- 2996 for the one with 20 NL 
links. In the middle pair where each one is at a temperature 
of 0.2 above its own critical temperature, the two g (r)’s 
are closer to each other and give rise to correlation lengths 
of the order of L, because the system is finite and is in the 
vicinity of phase transition. However in this case, one will 
see in the next section that there is some subtlety for the 
correlation function of the system with NL links. In the 
lower pair where each one is at a temperature of 0.4 below 
its own critical temperature, the two g (r)’s are almost the 
same and this is consistent with the fact that the correlation 
length should be finite at such temperatures. Therefore, 
even when the system is decorated with NL links, the num¬ 
ber of which is comparable to L, its correlation function is 
different but does not deviate very much from that of the 
regular lattice. 
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20 x 20 lattice with 20 non-local links 


20x20 lattice with 20 nonlocal links 




FIG. 4: comparison of Cy curves at 10 runs of simulation with all FIG. 5: Comparison of correlation functions at T = 2.3, which 

parameters the same except the configuration of non-local links. are obtained from 20 runs of simulation with same parameters, of 

a 20 x 20 Ising lattice with 20 non-local links. 


Effects of the Configuration of NL links 

So far we just add, at the beginning of each run of the 
simulation, a certain number of NL links on the lattice 
randomly based on the uniform distribution probability, 
without any further control of their locations in the lattice. 
Consequently, we may have different configurations of NL 
links on the lattice in different runs of the simulation with 
all other parameters fixed. Now, let us check if the con¬ 
figuration of NL links in the lattice has any impact on the 
system’s behavior. 

Actually, when the number of NL links is fixed, the Cy- 
curve is self-averaging. This is verified by simulation. As 
an illustration, Fig. 0 plots ten 6\/-curvcs of ten runs 
(arbitrarily taken from hundreds of runs) of the simulation 
with different NL link configuration but with all other pa¬ 
rameters the same. One can observe that all C\/-curves 
coincide very well and peak at the same T c . This suggests 
that for a system, Cy and hence T, are only sensitive to the 
total number of links (both regular and NL links) on the 
lattice. The reason is that Cy is a first order derivative of 
system’s energy with respect to temperature and the energy 
is only related the total number of links on the Ising lattice 
with nearest neighbor interaction. 

Nevertheless, if we look at the correlation function at 
temperatures not so far from the critical temperature, con¬ 
figuration of NL links does effect. Recall that we men¬ 
tioned the subtlety of the correlation function of a 20 X 20 
lattice with 20 NL links. This is illustrated in Fig. (0, 
where we compare correlation functions at T = 2.3 from 
20 runs (taken from hundreds of runs) of the simulation 
with different NL link configuration but with the same pa¬ 
rameters. The system is in the vicinity of the phase transi¬ 
tion, since T c ~ 2.5 in this case. Nineteen g (r)’s of them 
are similar in that they all indicate a correlation length of 
the order of L, with although different speed of falling off. 
We call them normal ones in the sense that they behave 
similar to the correlation function of the regular lattice at T 
of the same distance from the regular T c . However, one of 
them, namely the special one, falls off to zero rapidly and 


hence gives short correlation length of the system. This is 
indeed due to a special configuration of the NL links on the 
lattice. We show 20 g (r)’s in this figure, because the ra¬ 
tio of the special one to the normal one is roughly between 
1/30 and 1/20. 

Fig. © presents two different typical configurations of 
the 20 NL links on the 20 x 20 lattice. The configura¬ 
tion corresponding to the special correlation function in 
Fig. © is shown in Fig. 0i) (type A), where one can see 
that the whole upper left quarter of the lattice has virtually 
no NL links. On the other hand, the configuration in Fig. 
0) (type B), which corresponds to a typical one of the 19 
normal correlation functions, shows a much more evenly 
distribution of NL links. By doing many simulations, we 
found that the length l of NL links, measured by the num¬ 
ber of regular links taken from one end to the other end of 
the NL link, are statistically different in the two cases. For 
special ones, Z min > 5, and ( avg > 10, where ( min is the 
minimum length and / avg is the averaged l over the number 
of configurations of the same type. But in the normal cases, 
we have 2 < l min < 4 and ( avg < 10. Based on these ob¬ 
servations, this phenomenon can be explained as follows. 
Since we choose the uniform distribution probability of the 
NL links, in most of the cases one should obtain configu¬ 
rations of NL links of type B; this explains why there are 
more normal correlation functions than the special ones. In 
type B configurations, NL links are short and very evenly 
distributed, so their effect on the correlation function is 
well averaged out. As a consequence, the system is cor¬ 
related as there are no many NL links although T c is differ¬ 
ent. Whereas, in type A configurations, NL links are longer 
and packed in some region of the lattice. Provided this is 
so, the region where NL links gather are effectively a sub¬ 
system of much higher dimension with of course a much 
higher T' : , then the correlation function of this sub-system 
behaves as at a temperature much below T' c and gives rise 
to short correlation length as it should. Moreover, the ef¬ 
fect of NL links can hardly be averaged out in such cases. 
Therefore, one sees the special correlation functions. 
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FIG. 6: Two configurations of 20 non-local links on a 20 x 20 
lattice: (a) configuration corresponding to the special correlation 
function in fig. 0 ; (b) configuration corresponding to a normal 
correlation function in fig. © 


Conclusions and Future Works 

We conclude that 1) the Curie temperature of the system 
increases with the amount of NL links; 2) the Curie tem¬ 
perature is not sensitive to the configuration of NL links; 
3) normally, correlation functions of systems with NL links 
behave similarly to those of regular lattices, with noticeable 
difference in the decay rate; and 4) there exist some special 
cases where the correlation functions of systems with NL 
links behave very differently from regular ones due to some 
underlying special configurations of NL links on lattices. 
Although it is hard to exclude the boundary effect due to 
the small size of the systems we studied, the above results 
can still be valid for large systems, especially when we are 
only interested in the physics in some portion of the whole 
system. Nevertheless, to study large systems is of course 
demanding and is one of our next steps. Conclusions one 
and two suggest some manipulation of the interactions J 
in the model, because CV-curves and T c will be subject to 
the configurations of NL links if we choose a distance de¬ 
pendent interaction J (r). Conclusions three and four, on 
the other hand, suggest a modification on the distribution 
probability P (r) of the NL links. We see that correlation 


functions suffer from the configuration of NL links, so we 
may get a special type of correlation functions by control¬ 
ling the configuration of NL links through a well designed 
P (r). For example, we may take a J (r), which has been 
studied by others and has physical significance for a square 
spin lattice, and then look for a P (r) with the above prop¬ 
erty. This would imply that a physical system dominated 
by some interaction law can be equivalent to a system with 
non-local links, which however obeys a much simpler in¬ 
teraction law. Some of these may be done analytically. We 
can also extend the study to more complicated lattice sys¬ 
tems, or even to spin networks in the future. 
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